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We investigate the dynamics of micron-scale drops pushed across a hydrophobic or 
superhydrophobic surface. The velocity profile across the drop varies from quadratic 
to linear with increasing height, indicating a crossover from a sliding to a rolling 
motion. We identify a mesoscopic slip capillary number which depends only on the 
motion of the contact line and the shape of the drop, and show that the angular veloc- 
ity of the rolling increases with increasing viscosity. For drops on superhydrophobic 
surfaces we argue that a tank treading advance from post to post replaces the diffu- 
sive relaxation that allows the contact line to move on smooth surfaces. Hence drops 
move on superhydrophobic surfaces more quickly than on smooth geometries. 
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I. INTRODUCTION 

The question of how liquid drops move across a solid surface has long caught the interest 
of academic and industrial communities alike, with applications ranging from microfluidic 
devices to fuel cells and Inkjet printing. In many cases, efficient and effective control of 
the drop dynamics is highly desirable, and this relies upon understanding the internal fluid 
motion of the drop. 



It has been reported in the literature 







-|5j that a liquid drop may move in a variety 



of ways including sliding, rolling, tank treading and slipping. In some cases, a pearling 
instability may also be observed at the trailing edge of the drop [6]. Highly viscous drops 
roll rather than slide (sj , and on superhydrophobic surfaces drops appear to move very easily 



Using small partic 



es as tracers, it has recently become possible to access the velocity 



profile within a drop {2 5|. However it is not yet always clear how the internal fluid motion 



2 



is related to physical parameters such as fluid viscosity, the strength of the forcing, the 
equilibrium contact angle and hysteretic properties of the surface. Controlled experiments 
are typically possible only over a restricted range, and can be difficult for small drops of size 
below the capillary length. Moreover, certain parameters, such as surface roughness, can be 
difficult to control and yet may have an important role. Analytical calculations are possible. 
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11| contact angles. Computer 



but they are typically limited to small and high 

simulations are therefore highly desirable to bridge the knowledge from experiments and 
analytical theories. 

Several simulation methods have been developed to shed light on the dynamics of moving 



drops 
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12H181. The computational tool we use here is a 



20| solved using a lattice Boltzmann algorithm 2l|, |22| . This proves a useful approach 



mesoscopic diffuse interface model 



due to its ability to handle interfacial dynamics and complex geometries easily. For example. 



it has recent 



surfaces [24 



flows 
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y been used to study capillary fllling [23|, hysteresis of drops on patterned 



251], instabilities and detaching phenomena 15| and driven drops in external 



ITj . Diffuse interface models are mesoscopic approaches, which do not resolve the 



microscopic details of the contact line dynamics. As such, they provide a useful complement 
to molecular dynamics simulations. 

We concentrate on modelling micron-scale drops moving on hydrophobic and superhy- 
drophobic surfaces. We flnd that, for both hydrophobic and superhydrophobic substrates, 
the velocity proflle across the height of the drop falls into two regimes. Near the substrate 
the velocity varies quadratically with height corresponding to a sliding motion. This can 
be characterized in terms of a mesoscopic slip capillary number which depends on the equi- 
librium contact angle and the mobility of the contact line. Further from the surface the 
velocity varies linearly with height indicating rolling. There is no dependence of the angular 
velocity on microscopic details of the contact line motion, such as the mobility, but the 
rolling becomes faster with increasing viscosity for a given body force /. 

We identify several major differences between the motion on hydrophobic surfaces, which 
are smooth, and similar surfaces which are patterned with posts such that they become 
superhydrophobic. On the superhydrophobic substrates pinning of the contact line at the 
edges of the posts means that there is a threshold forcing below which the drop will not move. 
However, once they do start to move, they do so much more easily and drops can reach much 
larger velocities before being detached from the surface in agreement with the experimental 
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observations. We impose no-slip boundary conditions in the simulations, and on the smooth 
surfaces there is indeed no slip at the surface. However, the fluid velocity at the surface is 
non-zero and independent of the mobility for drops moving on superhydrophobic surfaces 
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28| suggesting that the interface moves by tank treading between successive posts, rather 
than by diffusion-mediated motion of the interface. 

The plan of the paper is as follows: In Sec. |TT]the numerical model is introduced, and in 
Sec. Illll the simulation set-up and the computational parameters are given. Considering first 
drops on smooth surfaces, in Sec. [IV]we describe the crossover from a Poiseuille-like flow to 
a pure rolling regime with increasing height and point out the relevance of the equilibrium 
contact angle. In Sees. |V]and|Vl]we investigate the dependence of the Poiseuille and rolling 
regimes on the viscosity and surface tension of the fluid and on the mobility of the contact 
line. In Sec. IVIII we compare the dynamics of drops on superhydrophobic surfaces. Finally, 
in Sec. IVIIIt we summarise and discuss our conclusions. 



We use a 



energy [IJI |20| 



II. THE MODEL 



binary fluid model, with components A and B, described by a Landau free 



^ = f{i^b + ^{dact>f)dV + ! ^sdS (1) 
Jv ^ Js 

where the bulk free energy density ipb is taken to have the form 

V'b = ynlnn + A(-i</.2 + l04). (2) 

The functional ([1]) is discretised on a cubic lattice with lattice spacing Ax, and At is the 
simulation time step, n is the total fluid density n = ha + ub, is the order parameter 
(f) = Ha — and c = Ax/ At. This choice of ipb gives binary phase separation into two 
phases with = ±1. 

The gradient term in Eq. ([1]) represents an energy contribution from variations in (p and 



is related to the surface tension between the two phases by 7 = a/SkA/Q and to the interface 



width through ^ = ^/k/A 20 |. 

The second integral in Eq. ([1]) is taken over the system's solid surface and is responsible 
for the wetting properties. The surface energy density is chosen to be ips = —heps [29], where 
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is the value of the order parameter at the surface. Minimisation of the free energy shows 



that in equihbrium the gradient in at the sohd boundary is 



20| 



-h. 



(3) 



The value of the phenomenological parameter h is related to the equilibrium contact angle 



9, 



eq V cq 



is taken with respect to the = 1 component) via 201]: 



h = V2kA sign (^-^cq) J cos (^^^ |l-cos 



(4) 



where a = cos ^ (sin^ ^cq) and the function sign returns the sign of its argument. 
The hydrodynamic equations for the binary fluid are 

dtn + da{nva) = , 
dt{nva) + dfj{nvavp) = -dpF^p + dii[7]{dpVa + daVii + dapd^v^)] + /„ , 



(5) 
(6) 
(7) 



Eq. ([5]) is the continuity equation, Eq. ([6]) the Navier-Stokes equation and Eq. ([71) the 
convection diffusion equation. M is a mobility coefficient. In Eq. ([6]) we have introduced 
the dynamic viscosity rj, and a bulk force for unit volume /„. The equilibrium properties of 
the fluid appear in the equations of motion through the pressure tensor and the chemical 
potential. Both can be obtained in the usual way from the free energy and are given by 



= A (-0 + 03) - kVV, 
"a/? = {Pb - 1(^70)^ - K,(f)d^^(f)) 6al3 + K{da(j)){dp(f)) 

p, = fn + A(-|02 + f04). 



(8) 
(9) 



To solve Eqs. ([5]), ([6]) and 



we use the free energy lattice Boltzmann scheme recently 
presented by Pooley et. al. [30|. The main advantage of this method is that the spurious 
velocities at the contact line are strongly reduced. This is particularly important when the 
two fluid components have different viscosities. 

A word of explanation is in order as to why we are using a binary fluid approach to 
describe a liquid-gas system. Recently several authors [23|, ISjJ, l32| have investigated capillary 
fllling using liquid-gas and binary lattice Boltzmann algorithms. They found that the binary 
model is able to reproduce the experimental results, but that the liquid-gas model produces 
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filling that is faster than expected. This is due to excessive condensation of the gas phase 
at the interface. The effect is suppressed only when the liquid density is much larger than 
the gas density, a regime difficult to access numerically. Therefore we use a binary model. 
However, our results arc equally applicable to a physical system where a liquid displaces 
a gas if the important physical parameters are the viscosities, rather than the densities of 
the fluid components. For convenience we will use liquid/gas terminology in the rest of the 
paper. 

The way in which the contact hue moves does depend on the choice of model. Here it 
is determined by the mobility M which determines how the diffusion across the interface 

is driven by variations of the chemical potential fi. Given that the model docs not include 
microscopic details we can say nothing about the mechanisms operating at the contact line 
in a physical system, and we regard M as a phenomenological parameter which controls the 
rate of slip. 

III. GEOMETRY AND SIMULATION PARAMETERS 

We consider drops of size smaller than the capillary length, so that we can ignore gravity, 
and sufficiently large that thermal fluctuations can be neglected. For a water drop this 
corresponds to a range of length scales between ~ 10~^m and 10^'^m. On these length 
scales the motion of the drop depends on its viscosity r] and surface tension 7 and a useful 
dimensionless measure of drop velocity v, allowing direct comparison to experiments, is the 
capillary number Ca — vrj/^. In the simulations we use three different values for the liquid 
viscosities {rj = 0.5, 0.833, 1.1667) and two for the surface tension (7 = 0.0267, 0.0533). The 
gas viscosity rigas is always chosen to keep the ratio rj/rjgas = 25 and a linear interpolation 
between 1] and r^gas is used to define the viscosity variation through the interface. 

A body force / is applied to the liquid phase. (For example, for a drop moving down an 
inclined plane with tilt angle a, f — n ■ g ■ sin a where g is the acceleration due to gravity.) 
/ is taken to vary with (1 + 0(r))/2 which provides a linear interpolation from / to when 
moving from the liquid to the gas. Other interpolating functions (e.g. tanh) were tested, 
without significant differences. Driving by the body force is normally reported in terms of a 
dimensionless Bond number Bo = f ■V/{'~f- Ly), where V is the volume of the drop, Ly is the 
length of the cylinder which is equal to the length of the simulation box in the y-direction. 
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FIG. 1: Simulation geometry for the patterned substrate. A cylinder of fluid is prepared in the 
suspended state and a force per unit volume f is applied along the x direction. 

and / is the body force defined in Eq. ([6]). We choose values of / between 0.05 ■ 10^^ and 
1.6 ■ 10-6 and 1/ = 11 ■ 10^ ■ Ly. 

Other parameters important in determining the drop motion are the mobihty M, where 
we use values from M = 0.125 to M = 2, and the equilibrium contact angle 6eq, considered 
below. 



We model cylindrical drops placed on a smooth surface or on a surface decorated by regular 
posts as shown in Fig. [TJ For the patterned substrates the use of cylindrical drops is a 
compromise between two dimensional simulations, which do not well describe the pinning of 
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the contact line on the posts (as they essentially replace the posts by ridges), and full three 
dimensional simulations which are prohibitively time consuming. 

For the smooth surfaces we mainly use a hydrophobic equilibrium contact angle 6eq = 
145°, but results at other contact angles are presented for comparison. For the patterned 
surface the posts were taken to have a square cross section of side 5* = 8 ■ Ax, height 
Dh = 10 ■ Ax, and to be equispaced hj = Dy = 9 lattice units in both plane directions 
(see Fig. [1]). A hydrophobic surface patterned by posts becomes super hydrophobic. There 
is a large increase in the contact angle 9c of a drop resting on the posts given by the Cassie 
formula cos 9c = <p cos ^cq + — 1 33[ where is the area fraction of the flat surface covered 
by posts. For the geometry we consider = S'^/{S + D^)^ ~ 0.22. Hence using 6'eq = 100° 
gives a Cassie angle 9c ~ 145° equal to the 6'eq used for the fiat surface geometry. 

As shown in Fig. [H the substrate is placed in the xy plane and the sliding motion is along 
the x-direction. Typically we used a simulation box of Lz = 180 X 1 X 180 

(lattice units) for the smooth geometry, and x Ly xL^ = 204 x 17 x 180 (lattice units) when 
posts were present. At small contact angles, for the fiat plane, was sometimes increased 
to Lx = 360 in order to minimize the interaction between the drop and its images created by 
the periodic boundaries conditions. Cylindrical drops were initialized with a radius R = 60 
centered 50 lattice units above the plane (or the top of the posts in the suspended states) 
and allowed to relax to equilibrium before the force was applied. 



IV. VELOCITY PROFILE INSIDE A DROP 



Understanding the dynamics of a moving drop is more difficult than for rigid objects 
because of the possibility of the fluid shearing. This means that, in general, the drop 
motion will include both sliding and rolling components. Fig. |2] reports a typical steady 
velocity profile found using the diffuse interface model Eq. ([T]) for an equilibrium contact 
angle ^eq = 145° in the laboratory and drop frames of reference. In the latter (Fig. [2]d) the 
presence of rolling is evident. 



In Fig. [3] we report the x-component of the drop velocity function of height 

z (computed along the line denoted H in Fig. [2]) for different applied forces /. There is 
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FIG. 2: (Color online) Velocity profile for a drop with contact angle ^eq = 145° (a) in the laboratory 
reference frame and (b) in the drop's center of mass reference frame. For clarity the velocities have 
been rescaled and we have not reported the flow of the surrounding gas. In panel (a) we illustrate 
the phenomenological description of the motion of a drop as superposition of a slip velocity Vgiip 
plus a rotational motion of angular velocity w. 

a clear separation into two regimes. For larger z the velocity scales linearly with height, 
corresponding to rolling dynamics. For smaller z the behaviour is close to quadratic and is 



well fit by a Poiseuille-like 3J] flow 



v^{z) = vo + —z -]--z^ , (10) 

7] 2 7] 

where H would correspond to the centre of the channel in the usual Poiseuille geometry and 
Vq allows for local slip at the liquid-solid interface. (For forces / < 0.5 • 10^^, where the 
drops are not too deformed, the quadratic term of the fit is in agreement with the z"^ term 
of Eq. (ITO!) to within 10% .) In our simulations we impose non-slip boundary conditions 
and therefore Vq = for smooth substrates by construction. We will see later that this is no 
longer the case for superhydrophobic surfaces. 

Experimental results are often fitted by simplifying the drop dynamics as a super- 
position of a rotation (with angular velocity w) and a constant sliding velocity VsUp- In this 
scheme the velocity profile in the x direction is a linear interpolation between VgUp and the 



9 



velocity at the top of the drop (see the arrows in Fig. [2^) 

Vx = Vshp + W- Z . (11) 




z [simulation units] 

FIG. 3: (Color online) Steady velocity in the x direction (f^) as a function of the height z. At 
large z a linear, rolling, regime is present (dashed) which crossovers a quadratic profile for small z 
(full lines). VsUp is defined as the extrapolation of the linear regime to z = 0. Data were obtained 
using r] = 0.833, 7 = 0.02667, 6'cq = 145° and M = 0.25. 



Guided by the experimental approach, we identify the angular velocity w as the slope of 
the linear part of the velocity profiles. Extrapolating the linear fit to z = (see the dotted 
curves in Fig. [3]) we define a slip velocity VgUp which characterizes the motion of the layers 
of fluid near the solid wall. (It is important to distinguish fsUp, and the microscopic slip 
velocity Vq which is the velocity at z = 0.) We will describe how VgUp and how the rolling 
region of the velocity profile, parameterized by u, depend on parameters such as viscosity, 
surface tension, mobility and the shape of the drop in Sections |V] and IVII respectively. 
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FIG. 4: (Color online) Variation of the velocity Vx with height z at different equilibrium contact 
angles but the same center of mass capillary number Ca = VcmV/l (where Vcm is the centre of mass 
velocity). Only for high equilibrium contact angles (^eq > 120°) is a linear region clearly present. 
The simulation parameters are the same as in Fig. [3j 



However first we present results showing that the equilibrium contact angle 6'eq strongly 
affects the crossover from the quadratic to the linear regime. In Fig. H] we report Vx for 
different equilibrium contact angles 70° < 6'eq < 150° at the same centre of mass capillary 
number Ca = where Vcm is the centre of mass velocity. For 6^,^ ^ 120°, a linear 

part of the plot cannot be distinguished; only for ^cq ~ 130° are the two regions clearly 
visible. The results are consistent with lubrication theory, valid for small contact angles, 
which predicts a quadratic profile jsS], and scaling calculations that predict rolling at high 



contact angle 
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11| . The applied force needed to produce the given capillary number in 



Fig. m increases with decreasing contact angle: for example the ratios are 1 : 1.7 : 2.7 for 
^eq = 70°, 105° and 150°. The dissipation is larger for drops of smaller contact angle because 
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FIG. 5: (Color online) Dimensionless slip velocity CusUp versus (a) the Bond number Bo and (b) 
cos 9pi — cos 9yi where Or and 6a are the receding and advancing contact angles. The curves depend 
on the mobility M and equilibrium contact angle ^eq but not, to with the precision of our data, 
the surface tension 7 or the viscosity r/. 



the region in which the velocity varies quadratically with height extends further from the 
substrate. 



V. SLIP VELOCITY, v,up 

We now discuss the dependence of the mesoscopic slip velocity on the parameters of the 
model. Fig. [5] shows a dimensionless measure of Vshp, CosUp = Vshp ■ rj/'j, plotted against 
a dimensionless measure of the forcing, the Bond number. Bo = f ■V/{'y ■ Ly). There is 
a data collapse onto curves which depend only on the mobility M and on the equilibrium 
contact angle ^eq, and Cosiip = /(^cq, M) ■ Bo holds well for Bo sufficiently small that the 
drop is not significantly deformed. As the mobility M increases fsUp increases for a given 
Bond number as expected, as it becomes easier for the contact line to relax. For smaller 
equilibrium contact angles VsUp decreases for a given Bond number as the dissipation of the 
drops increases. 
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Another property of the moving drop that is related to VgUp is the degree of deformation 
caused by the forcing. In Fig. [5]^b) we present data that suggest that, for small forcing, t>siip 
is linear in the uncompensated Young stress A = '-/{cosOji — cos 9 a) where Or and 6a are 
the receding and advancing contact angles. Again the slopes of the curves depend on the 
mobility and the equilibrium contact angles, but not on the viscosity or surface tension. We 
caution that the data here is noisy, as the deviations from the equilibrium contact angles 



3Q|. 



are small and the algorithm, for high viscosities, does not precisely reproduce ^eq 

It is plausible and consistent with our results that the mobility coefficient M affects CosHp 
only through the deformation of the droplets, in such a way that a more detailed scaling 
relation for CagUp could be written CosUp = a(6'eq) ■ Bo + (3{9eq,M) ■ A{Bo). Qian et al. 
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371] have also highlighted the importance of A in controlling interface motion. However 
the relation between our work and theirs is not yet clear as they consider a local microscopic 
slip near the interface, whereas our parameter Vgw-p follows from extrapolating from the bulk 
velocity profile. Work is in progress to investigate this further. 

VI. ROLLING 



The rolling motion displayed away from the surface by drops of higher equilibrium contact 
angle has a very different dependence on the model parameters. The rolling contribution 
to the velocity profile of the drop is w ■ z m. Eq. (ITT]) , corresponding to the linear part of 
the velocity profiles in Fig. [3l We estimate a typical rolling velocity as w ■ H, where H is 
a length comparable to the height of the droplet. A capillary number associated with the 
rolling motion can then be defined as Ca^^ = w ■ H ■ r]/'~f. For the graphs in Fig. |6]we have 
consistently used H = 100 lattice units. 

In Fig. |6]we plot Ca^, as a function of the Bond number for different mobilities, viscosities 
and surface tensions. Fig. [6^ shows that different viscosities do not collapse on a single 
master curve but rather that Ca^^ increases with rj. This means that the sliding and the 
rolling motion of the drops are affected in a different way by the viscosity. In particular the 
rolling component is favoured for high viscosity. This is not the case for the surface tension 
7: data for drops with different 7 (but equal rj) collapse onto the same curves to within the 
precision of the data. 
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FIG. 6: (Color online) Dimensionless angular velocity Ca^ = wqH/'j against (a) Bond number Bo 
(b) dimensionless slip velocity CosUp = rjVsHp/j. Curves for different surface tensions (but equal 
viscosity) collapse, while Ca^ increases with rj. In (b) the curves are also independent of mobility. 

In Fig. [6] Cayj is also larger for sets of data with high mobility. This is not unexpected 
because CosUp is also M dependent (Fig. [5^). To understand if the mobility affects the 
rolling part of the motion, we need to compare Ca^ at fixed CosUp. This is done in Fig. [6t^b) 
where it is apparent that different mobilities collapse onto the same curves showing that M 
affects the rolling part of the motion only through VsUp. 

For higher forcing the linear dependence between the capillary number and the Bond 
number is lost, co decreases, Vshp » w ■ H, and sliding dominates the motion. This 
corresponds to a regime where the drops are highly deformed and a further increase in Bond 



number causes them to detach from the substrate 15] . Similar behaviour has been observed 
in experiments for drops on hydrophobic surfaces. 

VII. DROPS ON PATTERNED SURFACES 



In this section we investigate the dynamics of a drop on a surface patterned with posts (see 
Sec. Illll and Fig. [T]for the geometric details). In particular we are interested in understanding 
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FIG. 7: (Color online) (a) The same as Fig. [3] but for superhydrophobic surfaces and for 0.29 < 
Bo < 0.54. The plot shows the velocity in the x direction (vx) as a function of the height z 
(measured from the top of the posts). At large z a linear, rolling, regime is present (dashed) 
with crossover to a quadratic profile for small z (full lines), (b) Vx for Bond number Bo = 0.371 
(/ = 0.9-10^^) and three different mobilities. Data were obtained using rj = 0.833 and 7 = 0.02667. 

analogies to, and differences from, motion across a smooth surface. Drops which lie on top 
of the posts in the suspended or Cassie state exhibit an apparent contact angle Oq given by 
Cassie's relation cos 6^0 = (p cos9eq + 4> — 1 |33|, where 6'eq is the equilibrium contact angle for 



the smooth plane. We consider cylindrical drops on surfaces patterned by obstacles which 
occupy an area fraction (p = 0.22 of the fiat surface, and a value ^eq = 100°, which gives 
6c ~ 145°, allowing a direct comparison with the results for a smooth surface in Sees. IVl 
andlVll 



In presence of posts the drop remains stationary for small Bond numbers 38l-l40l| . This 



occurs because the interface is pinned at the edges of the obstacles 4l|- For the geometry 
we consider here the depinning threshold Bot ~ 0.15. For higher Bond number we observe 
oscillations in the velocity which reflect the position of the drop relative to the posts and 
become smaller further from the pinning threshold Bot- The steady velocities were therefore 
computed by averaging over a time interval much longer than the time required for the drop 
to cross a post. 
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In Fig. El^a) we plot the velocity profile as a function of the distance from the top of the 
posts for Bond numbers 0.29 < Bo < 0.54. This plot should be compared to Fig. |3] for 
the smooth surface. Just as for the smooth surface, the velocity profile exhibits a crossover 
from a quadratic to a linear regime. The main difference is the presence of a microscopic 
slip velocity at 2; = 0, defined as Vq in Eq. ( [T0|) . for the superhydrophobic substrate [27|, |28 |. 
These results show that a crossover from a quadratic to a linear dependence of velocity on z 
is not dependent on having zero Vq but is also present if some slip occurs at the solid-liquid 
interface. We remark that it is difficult to obtain a precise value for f from the simulations 
because of the presence of spurious velocities at the surface. 

An important difference between the smooth and the patterned geometry is illustrated 
in Fig. El^b) where we plot against z for Bo = 0.371 and different values of the mobility 
M. The velocity is unaffected by the mobility, in stark contrast to the smooth surface where 
Cttsiip depends strongly on M (compare Figs. [5]). This provides evidence that the contact 
line is moving onto successive posts using a tank-treading mechanism, rather than through 
relaxing the interface distortion by diffusion. 



Fig. [8] is a plot comparing the centre of mass capillary number as a function of Bond num- 
ber for drops on hydrophobic and superhydrophobic surfaces. Full lines refer to the smooth 
surface, while the broken line is for the patterned geometry. As before, the equilibrium 
contact angle of the smooth surface is ^^cq = 145°, while for the patterned plane ^cq = 100° 
corresponding to a Cassie angle of 6c ~ 145°. For Bo > 0.3 the drop on the smooth plane 
is very deformed from its equilibrium shape and a small further increase in Bo leads to the 
drop detaching from the substrate. For the drop on the posts steady motion is only possible 
above a threshold Bond number due to pinning, but then the drop moves faster and with 
much less deformation. Thus it is possible to push the drop much harder and to achieve 
velocities larger by a factor ~ 4 before detachment occurs. 

Increasing the bulk force / leads to a drop that detaches from the surface. We have 
observed that, for superhydrophobic surfaces, there is a long transient regime in which the 
droplet accelerates before it finally detaches. This was not observed in the case of smooth 
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FIG. 8: (color online) (a) Drop center of mass capillary number Cflcm versus Bond number Bo 
for drops on hydrophobic and superhydrophobic surfaces. The equilibrium contact angle of the 
smooth surface is 6eq = 145°, while for the patterned plane it is 6eq = 100° giving a Cassie angle 
9c ~ 145°. (b) Steady state shapes of the drops at different Bo. 
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planes. Similar accelerating drops have been observed in experiments on superliydropliobic 



surfaces 



3| although the acceleration measured in the simulations is much smaller (0.2% x /) 
than that in the experiments (~ /). 

VIII. DISCUSSION 

We have used a mesoscopic simulation approach to investigate the velocity profile of 
hydrophobic and superhydrophobic drops subject to a constant body force. 

For both hydrophobic and superhydrophobic surfaces we found that the velocity profile 
perpendicular to the substrate Vx{z) comprises two regions. Close to the surface Vx{z) is 
quadratic in z, as in a Poiseuille flow. This is the velocity profile found using lubrication 



theory 



35| . Further from the substrate Vx{z) is linear in z describing a rolling motion, as 



predicted for drops of large 6'eq by scaling arguments 10, lUl ■ The crossover is clearly seen for 
equilibrium contact angles ^ 120°, for smaller contact angles (e.g. 70°) the velocity profile 
is quadratic for all z. 

Fitting the linear portion of the curve as 

Vx = Vslip + w ■ z (12) 

defines a mesoscopic slip velocity WsUp which characterizes the quadratic region of the flow 
field. We find that CosHp is independent of the viscosity and the surface tension but depends 
on the the shape of the drop, through the equilibrium contact angle, and on the mobility 
parameter M which controls the ability of the contact line to move across the surface. This is 



reminiscent of the 'inner' region, identified by Cox 42|, which encapsulates the microscopic 
physics. 

The rolling part of the flow profile was found to be independent of the mobility. Thus 
it might be identified as the 'outer' region of the flow which is, as expected, independent of 
the microscopic details of the contact line motion. The angular velocity is independent of 
surface tension, as expected, but increases with increasing viscosity. 

A similar behaviour is seen for drops on superhydrophobic surfaces, but there is also a 
microscopic slip velocity, that is, a non-zero velocity of the fluid adjacent to the surface. We 
stress that this occurs despite the no-slip boundary conditions imposed in the simulations. 
It suggests that the contact line moves not by diffusion, but by tank-treading from one post 



18 



to the next. Further evidence for this is that, for the superhydrophobic surfaces, VsUp is 
independent of the mobihty. 

Pinning on the posts leads to contact angle hysteresis and it takes a non-zero force to 
initiate drop motion. Once the drop is moving it can be pushed more easily than on a 
flat surface and can reach higher velocities before flying off the substrate. The drop slowly 
accelerates before leaving the surface. 

Rolling dynamics have been observed in experiments. In jsl an almost pure rolling motion 
was measured for glycerol drops with a superhydrophobic coating (6'eq = 165° with an 
hysteresis A^^ < 10°) and diameters between 1.2 and 8 mm. The dominance of rolling 
was ascribed to the high viscosity of glycerol (950 times bigger than water) and the high 
equilibrium contact angle. Suzuki et al. reported a linear profile for the velocity inside 
water drops with 6'eq ~ 110° but did not resolve any quadratic regime. However slip at the 
solid substrate was observed, and identified as fsUp, and three different cases were reported 
in which f sUnWas negligible, comparable and dominant with respect to the rolling part of the 
motion. In 3| predominantly sliding drops were observed in the superhydrophobic regime 
(^cq = 150°), while for lower 6'eq rolling motions were also resolved. 

A recent molecular dynamic simulation jl^ gave a Poiseuille profile ffTOj) over all the 
height of the drop. In particular no rolling regime was observed. Possible explanations 
for the difference between these results and ours are that the equilibrium contact angle in 
the molecular dynamics simulations was 130°, corresponding to a region in our simulations 
where rolling was not well resolved, that molecular dynamics typically accesses much smaller 
length scales than the diffuse interface approach or that the intrinsic slip mechanism at the 
contact line was different in the two methods. 

Clearly more work needs to be done to unify experiments, mesoscopic simulations and 
microscopic simulations to give a clear picture of how the way in which a drop moves across 
a surface depends on its geometry, the parameters describing the fluid and the way in which 
the drop interacts with the substrate. To work in this direction our next aim is to gain 
a better understanding of how the size of the drop and the degree of intrinsic slip at the 
surface affect the simulation results. 
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